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Abstract 

We investigate the asymptotic behavior of the eigenvalues of the sum 
An + U^BnUn , where An and Bn are deterministic N x N Hermitian 
matrices having respective hmiting compactly supported distributions ^ 
and v, and Un is a random N x N unitary matrix distributed according to 
Haar measure. We assume that An has a fixed number of fixed eigenvalues 
(spikes) outside the support of /i whereas the distances between the other 
eigenvalues of j4jv and the support of fi, and between the eigenvalues 
of Bn and the support of v uniformly go to zero as A'^ goes to infinity. 
We establish that only a particular subset of the spikes will generate some 
eigenvalues of j4jv + UnBnUn outside the support of the limiting spectral 
measure, called outliers. This phenomenon is fully described in terms of 
free probability involving the subordination function related to the free 
additive convolution of /x and u. Only finite rank perturbations had been 
considered up to now. 

1 Introduction 

The set of possible spectra for the sum of two deterministic matrices Ajsf and 
Bn depends in complicated ways on the spectra of An and Bn (see [H]). 
Nevertheless, if one adds some randomness to the eigenspaces and assumes 
them to be in generic position with respect to each other, when N becomes 
large, free probability provides a good understanding of the global behavior of 
the spectrum of the sum of matrices. Indeed, ifXjv — An + U^BnUn, where 
C/at is a Haar unitary random matrix (i.e. from the set of unitary matrices 
equipped with the normalized Haar measure as probability measure), if the 
spectral measures of An and Bn converge weakly towards respective compactly 
supported distributions fj, and then building on the groundbreaking result of 
Voiculescu , Speicher proved in [35] the almost sure weak convergence of the 
spectral measure of Xn to the free convolution /i ffl z^, which is known to be a 
compactly supported probability measure on R. We refer the reader to [55] for 
an introduction to free probability theory. 

In the authors investigated the case where An has a finite rank r 

independent of N. More precisely, they considered a deterministic Hermitian 
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perturbation matrix An having r non-zero eigenvalues 71 > • • • > 7s > > 
7s+i > • • • > 7r- Note that in that case, fi = So and the global limiting 
behavior of the spectrum oi is not affected by such a matrix A^. Thus, 
the spectral measure of Xn = An + U^BnUn still converges to the limiting 
spectral measure of Bn- Nevertheless, in [TT], the authors uncovered a phase 
transition phenomenon whereby the limiting value of the extreme eigenvalues 
of Xn differs from that of Bn if and only if the eigenvalues of An are above a 
certain critical threshold: 

Theorem 1.1 (Theorem 2.1 in [TT]). Denote by XiiXN) > ••• > A7v(^Ar) 
the ordered eigenvalues of Xn- Let a and b be respectively the infimum and 
supremum of the support of v. Assume that the smallest and largest eigenvalue 
of Bn converge almost surely to a and b. Then, we have for each 1 < i < s, 
almost surely, 



while for each fixed i > s, almost surely, Xi(XN) —^n^+oo b. Similarly, for the 
smallest eigenvalues, we have for each < j < r — s, almost surely, 



while for each fixed j > r — s, almost surely, XN-j{XN) ^n^+oo a. Here, 



is the Cauchy-Stieltjes transform of v, Gjj ^ is its functional inverse. 

Note that |TT] lies in the lineage of recent works studying the influence of 
some finite rank additive or multiplicative perturbations on the extremal eigen- 
values of classical random matrix models, the seminal paper being [5] where 
Baik, Ben Arous and Peche pointed out the so-called BBP phase transition (cf. 
[25l[5l[6] for sample covariance matrices, [22l |28l [HI [151 [29] for deformed Wigner 
models and [2 7) for information-plus- noise models). Such problems were first 
extended to non-finite rank perturbations in [30] and [3] for sample covariance 
matrices and in [16] for deformed Wigner models. In this last paper, the authors 
pointed out that the subordination function relative to the free additive con- 
volution of a semicircular distribution with the limiting spectral distribution of 
the perturbation plays an important role in the fact that some eigenvalues of the 
deformed Wigner model separate from the bulk. Note that in ^14j, the author 
explained how the results of [3D] and [5] in the sample covariance matrix setting 
can also be described in terms of free probability involving the subordination 
function related to the free multiplicative convolution of a Marchenko-Pastur 
distribution with the limiting spectral distribution of the multiplicative pertur- 
bation. 
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In this paper, we investigate the asymptotic spectrum of the model Xjv = 
An + U^BnUn without the extra requirement that one of the limiting spec- 
tral measures fj, and v of the respective deterministic Hermitian matrices 
and Bn is a point mass. We assume that An has a fixed number r of fixed 
eigenvalues (spikes) outside the support of /j. whereas the distances between the 
other eigenvalues oi An and the support of n, and between the eigenvalues of 
Bn and the support of v uniformly go to zero as N goes to infinity. We are 
interested in the following questions: are some of the eigenvalues of Xn almost 
surely separating from the bulk, that is being located outside the support of the 
limiting spectral distribution nSvl What is the characterization of the spiked 
eigenvalues of An that generate such outliers in the spectrum of X^v? Is there 
an interpretation of these phenomena in terms of the subordination functions 
related to the free additive convolution of /i and 

In the particular case where r = 0, it is proved in [15] that, almost surely, 
for large enough iV S N, all eigenvalues of Xn are included in a neighborhood of 
the support of the limiting spectral distribution fiWv. In the following, we will 
therefore assume that r > 1. Thus, this paper may be seen as an extension of 
since it extends the framework of to non-finite rank perturbations but 
also as an extension of jl6| since it extends the free probabilistic interpretation 
of outliers phenomena in terms of subordination functions described in 1161 for 
Wigner deformed models to deformed unitarily invariant models. Here, the 
characterization in terms of subordination functions of the spiked eigenvalues 
of An that generate outliers in the spectrum of Xn turns out to be more 
complex that the one presented in [16] for deformed Wigner models, but it is a 
completely natural extension as we explain in Remark 1 3. 2 1 It is worth noticing 
that we uncover here a new phenomenon: a single spiked eigenvalue of An may 
generate asymptotically a finite or countably infinite set of outliers oi Xn- This 
comes from the fact that the restriction to the real line of some subordination 
functions may be many-to-one, unlike the subordination function related to free 
convolution with a semicircular distribution studied in [TB] . 

The approach of the proof of our main result (i.e Theorem I4.ip is in the 
spirit of [11] and comes down to prove the almost sure convergence of a certain 
r X r matrix, involving the resolvent of the deformation of W^BnUn by some 
matrix A'j^ without spikes. Its almost sure convergence is proved by establish- 
ing an approximate matricial subordination equation and using a concentration 
argument. Then, the problem is reduced to solving an equation involving the 
spikes and the subordination function related to the free convolution of /i and 

V. 

The paper is organized as follows. In Section 2, we introduce the additive de- 
formed models we consider in this paper; we also introduce some basic notations 
that will be used throughout the paper. Section 3 is devoted to definitions and 
results concerning free convolution and subordination functions, some of them 
being necessary to state our main result Theorem 14.11 in Section 4. The proof 
of Theorem 14. II is presented in Sections 4 and 5. More precisely, we explain in 
Section 4 how the proof comes down to prove the almost sure convergence of a 
r X r matrix and Section 5 deals with the proof of this convergence. 
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2 Notations and presentation of the model 



Throughout this paper, we will use the following notations. 

- C+ will denote the complex upper half-plane {z G C, > 0}. Similarly, 
C" will stand for {z eC,Qz < 0}. 

- Wc will denote by M„i(C) the set of m x m matrices with complex entries 
and GLm(C) the subset of invertible ones. || || will denote the operator 
norm. 

- For any matrix M, we will denote its kernel by Ker(M). 

- Eij stands for the matrix such that {Eij)ki — SikSji. 

- For any N x N Hermitian matrix M, we will denote by 

Ai(M) > ... > Ajv(M) 

its ordered eigenvalues. 

- For a probability measure r on M, we denote by supp(r) its topological 
support. 

- C,Ci,C2,C' denote nonnegative constants which may vary from line to 
line. 

In this note, we consider the following model Xjv = ^jv + U^B^Un, where: 

• An is a deterministic N x N Hermitian matrix whose spectral measure 
/iAiv •= jf^iLi^XiiAn) weakly converges to some compactly supported 
probability measure /j. on M. We assume that there exists a fixed integer 
r > (independent from N) such that has N — r eigenvalues a^-^' 
satisfying 

max dist(Q;*;-^'', supp(/i)) — > 0. 

We also assume that there are J fixed real numbers 0i > . . . > 9j inde- 
pendent of N which are outside the support of /x and such that each 9j is 
an eigenvalue of An with a fixed multiplicity kj (with J2j=i — '')• The 
9j's will be called the spikes or the spiked eigenvalues of An- 

• Bn is a deterministic N x N Hermitian matrix whose spectral measure 
A*Bjv •= ;^ Sil=i '^Ai(Bjv) weakly converges to some compactly supported 
probability measure u on M. We assume that the eigenvalues jSj^^ of Bn 
satisfy 

max dist(/5!^^\ supp(i/)) — > 0. 

l<j<N JV^-oo 

• Un is a random N x N unitary matrix distributed according to Haar 
measure. 
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3 Free convolution 



Free convolution appears as a natural analogue of the classical convolution in 
the context of free probability theory. Denote by A4 the set of Borel probability 
measures supported on the real line. For fj, and v in A4 one defines the free 
additive convolution ffl z/ of /i and i' as the distribution oi X + Y where X and 
Y are free self-adjoint random variables with distribution fi and v. We refer the 
reader to [31] for an introduction to free probability theory and to |33] and [TU] 
for free convolution. In this section, we recall the analytic approach developed 
in |33j to calculate the free convolution of measures, we present the important 
subordination property and some related fundamental results we will refer to 
later. 



3.1 Additive Free convolution 

For any positive finite Borel measure r on M, the Cauchy-Stieltjes transform of 

T 

: C \ supp(t) ^ C, G^(z) = / 

is analytic, maps the upper half-plane C+ into the lower half-plane C~ and satis- 
fies the conditions Gr{'z) — Gt{z), hm^^^+oo iyGr{iy) = t(E). These conditions 
in fact characterize functions which are Cauchy-Stieltjes transforms of positive 
finite measures. For probability measures t with compact support, Gr is ana- 
lytic on the neighbourhood of infinity {z G C: \z\ > max{|a;|: x £ supp(r)}}. 
The Cauchy-Stieltjes transform allows us to recover the measure r as the weak*- 
limit ^ 

dr(x) = lim — 3G,-(a; + iy). 

For the absolutely continuous part of r, the situation is better: the relation 

4^ = lim — -"^sGrix + iy) 
ax y^o TT 

holds Lebesgue-a.e. as equality of functions. For some purposes, it is convenient 
to work with the reciprocal Cauchy-Stieltjes transform, which is the analytic 
self-map of the upper half-plane defined by: 

The Cauchy-Stieltjes transform of a compactly supported probability measure 
r is invertible in the neighborhood of infinity, with functional inverse, denoted 
by G^^, defined in a neighborhood of 0. Define then the R-transform of r by: 

Rr{z) = G:;\z)--. 

z 
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Given two compactly supported probability measures fi and i^, there exists a 
unique probability measure A such that 

^^A = Rn + Rv 

on a domain where these functions are defined. The probability measure A is 
called the free additive convolution of /i and v and denoted by /i ffl i/. The 
support of the probability measure fi'Siv being a compact set, we will denote it 

by 

K :— supp(/i ffl v). 
Moreover, given e > 0, we will use the following notations: 

:= {a; e M I d{x, suppin S i^) < e}, 

iCf := {z e C I d{z, supp(Ai ffl i^) < e}. 

3.2 Free subordination phenomenon 

We recall the subordination phenomenon for the Cauchy-Stieltjes transform of 
the free additive convolution of measures. Given Borel probability measures fj, 
and on M, the Cauchy-Stieltjes transform of the free additive convolution /xfflz/ 
is subordinated to the Cauchy-Stieltjes transform of any of /i or v. there are 
two analytic self-maps of the upper half-plane wi,W2 : C"*" C+ such that: 

Vz e C+, G^ffl,(z) = G^{uji{z)) = G.{uj2{z)). 

The subordination maps wi,aj2 are also related by: 

VzeC+, uji{z)+uj2{z)^z + F^siAz)- (3.1) 
Moreover, for j e {1,2}, 

hm ^ = 1. 

y— H-oo iy 

It then follows, using Nevanlinna representation of analytic self-maps of the 
upper half-plane, that: 

equality can occur only when one of the measures iJ,,v is a point mass (more 
specifically, 3z e C+, 9a;i(2:) = 9z Vz € C"*", 3a;i(z) = 3z ly is a 

point mass, and 3z £ C+, 3a;2(z) — 5z ■^=> Vz e C^, 3w2(z) = 3z <^=» ^ 
is a point mass). These results, first obtained in full generality in |13) , have been 
given a new interpretation in terms of Denjoy- Wolff points of analytic functions 
in [7] . We will state the upper half-plane version of the Denjoy- Wolff theorem 
below: 

If / : C+ — > C+ is analytic and not a Mobius transformation, then only one 
of the following three cases can occur: 
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1. there exists a unique point cj G C"*" so that /(w) — uj and |/'(a;)| < 1. 
Then the iterations /°" converge uniformly on compacts to the constant 
function taking the value uj; 

2. there exists a unique a; G M so that linv^o f{ii~ + co) — uj and 

f(ir + Lu) — id 
< lim ^ — < 1. 

rlQ ir 

Then the iterations /°" converge uniformly on compacts to the constant 
function taking the value uj; 

3. limr-f+oo — oo and 

1 < lin, ZM < oo. 

rt+oo ir 

Then the iterations /°" converge uniformly on compacts to infinity. 

The point uj in situations 1. and 2. (and infinity in 3.) is called the Denjoy- Wolff 
point of /. Note that the second limit in 2. above always exists in (0; +oo] as 
soon as lim^^o /(*^ + uj) = uj, and is called the Julia-Caratheodory derivative of 
/ at UJ. Denjoy and Wolff proved that any analytic function / : C+ — > C+ as 
above has a Denjoy- Wolff point (see ^\ for details). Using this result, a new 
proof of Biane's subordination result for free convolution was given in [7], by 
identifying uJi{z) as the Denjoy- Wolff point of the function 

f,{uj) F^{Ff,{uj) ~ UJ + z) - {F^{uj) - UJ + z) + z. 

We will need the following lemma collecting results on extensions of the 
subordination maps: 

Lemma 3.1. For j G {1,2}, the function ujj, defined on C^, has an extension 
(still denoted by ujj) to C so that: 

(a) ujj is continuous on C+ U M; 

(b) a;i({oo}UM\supp(^fflt/)) C {cx)} U R \ supp(/i); 
(b') a;2({oo}UM\supp(^ffli/)) C {oo} U R \ supp(^); 

(c) Vz e C \ R,ujj{z) ^ ujj{z); 

(d) ujj is meromorphic on <C\ supp(/x ffl v). 

Proof. As noted in [S] Theorem 3.3], ujj\f^+ has a continuous extension to the 
real line. This proves (a). Letting z tend to x e M \ supp(/i ffl ly) in (|3.ip . 
one notices that ujj{x) necessarily belongs to the boundary of C"*", because so 
does F^^^{x). To be more specific, if F^Hy(a;) G R, then ujj{x) G K by relation 
(|3.ip and the Nevanlinna representation. If F^gj^{x) = oo, then x is an isolated 
simple pole of F^si^, and the same (|3.ip tells us that cither uji{x) = oo and 
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i^2{x) = — mi(/i) + X, or ^2(2:) = 00 and uii{x) = —mi{v) + x, and thus x is 
an isolated simple pole of exactly one of the two ujj. (Here mi(T) denotes the 
first moment of r.) This proves (b), (b') and (d). We then extend Uj to C~ by 
requiring: 

Vz e ,LUj{z) = Li>j{'z), 

as in [13], so that (c) holds. □ 

It should be noted that the proof from [7] shows that uji{z) is the Denjoy- Wolff 
point of fz only for z S C'^. Indeed, it is still an open problem whether this 
is true for all z E C+ U K. However, if a: S M is so that wi is analytic in x 
and real on some interval around x, the Julia-Caratheodory Theorem requires 
u!'i{x) G (0,+c»). As noted in the above lemma, this implies that is real 
meromorphic around uJi{x) and real analytic around W2(a;) — -F)i(wi(a;)) — 
u}i{x) + X. Thus, fx will be real analytic around uJi{x). Taking limits shows 
that fxiuJi{x)) = uJi{x) and uj[{x) = dxfxi^i{x)) [1 - du^fxii^iix))]'^ ■ Since 
dxfx{i^i{x)) = Fl{u!2{x)) > 0, in order for uj[{x) > as required by the Julia- 
Caratheodory Theorem, we must have duifx(.<^i(x)) < 1, and hence w ^ 
has ijJi{x) as a fixed point in which the derivative is less than one. This is 
the case 2. of the upper half-plane version of the Denjoy- Wolff theorem given 
above. Thus, in this case uji{x) is still necessarily the unique Denjoy- Wolff point 

of W 

Let us give a slightly different formulation for the results of 7\ concerning the 
subordination functions, more appropriate to the needs of our paper. Assume 
neither nor 1/ is a point mass and denote 

= - z, K{z) = Fy{z) - z. 

We re- write the ideas of 0, where wi and 0J2 are identified as the Denjoy- Wolff 
points of self- maps of C"*" indexed by z (see above), but with a formulation 
chosen to make the statement somehow independent of complex dynamics. It 
follows from [7] that wi and uj2 are identified uniquely by the following system 
of equations 

r CJi(z) - hy{uj2{z)) = z 
\ uj2{z) - h^{uji{z)) = z 

We can look upon this system as an implicit equation for a two- variable map: 
we define /(wi, W2, z) = (wi — h^(w2) — z,W2 — /i^t(it;i) — z). A straightforward 
application of the implicit function theorem indicates that h'^{wi)h'^{w2) = 1 
is the only obstacle to an analytic solution z {LUi{z),i02{z)) to the equation 
f{uJi{z),uj2{z), z) — (0,0). The analysis in [7 shows that 

VzeC+, \h'^{u;,{z))hl{u;2{zj)\ < 1. 

This is a consequence of the fact that non-trivial - i.e. not Mobius - maps 
have derivatives less than one in their Denjoy- Wolff points. Since the func- 
tions h and uj map the parts of R which lay in their domains of analytic- 
ity in R and their derivatives on R are necessarily positive, we conclude that 
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< h'^{ijJi{x))h'i,{uj2{x)) < 1 for any a; G R in the domain of analyticity of uji and 
0J2- In particular, < h'^{uji{x))hl{hf^{uji{x)) + a;) < 1, and the obstacle to an 
analytic solution wi is described by the equality h'^{u!i{x))h'^{h^{u!i{x))+x) = 1. 

Theorem 3.1. Assume that neither /i nor v is a point mass. Given e M\ 
supp(/i), then p G M is a solution of the equation 

uji{p) = 9 

belonging M \ supp(^ ffl v) if and only if h^{9) + p belongs to the domain of 
analyticity of h^, and p is a solution of: 

h^{h^i0) + p) - 9 + p = 0, < h'^{0)hl{h^{0) + p) < I. (3.2) 

Proof: Let p E M. \ supp(/i ffl i^) be a solution of the equation uJi{p) — 0. 
Taking the limit z — > p in the second equality of the system above, we obtain 

hf_,{9) + p = oj2{p) 

which belongs to {oo} U M \ supp(i/) by Lemma [??T] (b'), and therefore to the 
domain of analyticity of /i^. Taking the limit z — )• p in p.ip . one gets 

h„{h^{9) + p) = h^{uj2{p)) = uJi{p) - p = 9 - p. 

Since p € R\supp(/ifflj^), uji is analytic at p (see Lemma [XT] (d) - by hypothesis 
wi(p) = 9 ^ oo), and there must be no obstacle to the existence of an analytic 
solution around p as in the discussion above. Hence 

0<h'^{9)h'M9) + P)<^- 

Conversely, let p G M be such that h^{9) + p belongs to the domain of ana- 
lyticity of hi^, and p is a solution of (j3.2l) . This in particular implies that 9 is 
the Denjoy-Wolff point of fp{uj) = /i^(/i^(a;) + p) + p, w e C+, as fp{9) = 9 
and fp{9) g]0,1[. An application of the implicit function theorem shows that 
the dependence of 9 on the complex variable p is analytic around the given real 
point, and has a positive derivative. Analytic continuation, the uniqueness of 
the Denjoy-Wolff point and the considerations above guarantee that uJi{p) = 9. 
□ 

Remark 3.1. We should note that the set of points p that satisfy the equality 
from p.2p . being the set of zeroes of an analytic map, is necessarily discrete, by 
the principle of isolated zeroes. Therefore the set of solutions of uji (p) — 9 in 
K. \ supp(/i ffl v) is discrete as well. This set of solutions may be empty, finite 
or countably infinite, as illustrated in the next remarks. 

Remark 3.2. If v is ^-infinitely divisible, it is known (see 1121) that lui is a 
conformal bijection between C+ and a simply connected domain 51 C whose 
inverse is the restriction to 51 of the map H defined on C+ by: 

H{z) -.^ z + R,{Gp{z)), 
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where R^, is the R-tranform of the measure v. Under this extra hypothesis, the 
content of Theorem \3.1\ reads simply: p € R is a solution of the equation 

Loi{p) = e 

belonging io R \ supp(/i ffl v) if and only if 

P = H{9), H'{9)>0. (3.3) 
Indeed, i/ p e R satisfy p.3p , notice first that 

h^{e) + p = F^{e) + R,{G^{d)) = G-\G^{e)), 

which is in the domain of analyticity of hi, (recall that 6 ^ supp(/i) by hypothesis 
and thus Gu{h^[9) + p) = G^{9)). Then 

h..{h^{e) +p)+p = F„{G-\G^{e))) - h^{9) = F^{9) ~ h^{9) = 9. 

Finally, because of the relation 

H{9) = 9~h,{h^{9)+p), 

the condition H' [9) > implies the inequality in (13.21) . 

Conversely, assume that p G R is such that h^{9) + p belongs to the domain of 
analyticity of hy and satisfies p.2p . The condition p.3|) being equivalent to 

9 = p-R,{G^bAp)), H\9)>0, (3.4) 

it is possible to check that p— Ru{G^^,^{p)) is the unique D enjoy- Wolff point 
of fp, which yields the equality above. The inequality follows immediately from 
the inequality in p. 21) . Note that there is at most one solution belonging to 
R \ supp(/i ffl v) which satisfies the equation 

L0l{p)^9 

when V is ^-infinitely divisible. 

Example 3.1. We should note that there are many examples in which the re- 
striction to the real line of subordination functions are many-to-one on their 
domain of analyticity. A trivial example comes from free Brownian motions, 
in which the spikes are associated to the matrix approximating the semicircular 
distribution (note the difference from the problem studied in U^f !). Consider a 
Bernoulli distribution b = ((5_i + <5i)/2 and a standard semicircular 7f of vari- 
ance t. We know that Gi,fs-yt{z) — G-y^{uji{z)) = Gh{^2{z)), where lj2 satisfies 

uj2{z) ^ z -tGb(uj2{z)), zeC+UR. 

We claim that for t > Q small enough, there are real values taken twice by uji . 
Indeed, 

— 1 —2 
hbih.y, (9) + p) = - — T-TT = , , 

^ "^'^ ^ h-,M+P -9+s/W^t + 2p 
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and then 



hb{h^,{6)+P) =S-P ^ + (^612 -4t~ W)p + 0^- 6^/6^ - At 
The two solutions are 



361 - Vd^ - At ± \/2e^ + 29V0^ -At + 16 
Ph2 ^ . 

For these solutions to be real, we must have \6\ > and 



e'^+8-2t> -Oy/e^ - At. 

For our example we shall consider a positive spike, and so the requirements 
translate into 9 > 2\Jt and 6^ > 2t — 8. For t < 1, this reduces to 9 > 2\ft. On 
the other hand, the condition h'^^{9)h'^{h^f^{9) + p) < 1 means 

KM 

We recall from above that -g — j^y+p = ^ ^ P- Thus this inequality is simply 

h'^^{0){p - 9f < 1. We note that h'^^{0) = ^^73= - 5 G]0, 1[ for all > ZVt. 

Moreover, linig_>._|_oo — and, again when — > +00, the two values of 

p tend to infinity at a speed of the order of 9 and at zero with a speed of order 
^ , respectively. Now we conclude easily: for 9 large enough and t < ^ ( strict 
inequality!), both conditions in p.2p are satisfied for the existence of p, and thus 
we obtain two spikes explicitely computed. 

Remark 3.3. In the case, not covered by Theorem \3.1\ where p, or v is a point 
mass, UJi is equal either to F^, or to a real translation. In the second case, the 
equation 

LJl{p)=9 

is trivial. In the first case, Fy establishing continuously increasing bisections 
between each connected component of the complement of the support of v and 
some subsets of R, the set of solutions may be determined by real analysis: 
for example, if v — ^2~"(5i/„, then there are countably many open intervals 
included in ]0, 1[ which are mapped bijectively onto M by F^. 

In this paper we shall be concerned almost exclusively with wi, so for sim- 
plicity we shall adopt the 

Notation convention: ui = ui. 

4 Main results and sketch of proof 
4.1 Main result and examples 

Definition 4.1. For each j G {1,..., J}, define Oj the set of solutions in 
R \ supp(/i ffl v) of the equation 

uj{p) = 0,, (4.1) 
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and 

U o,. 

i<j<J 

Recall that the sets Oj defined above may be empty, finite, or countably infinite. 

Theorem 4.1. Denote by sp{X) the spectrum of the operator X . The following 
results hold almost surely: 

• for each p € Oj, for all small enough e > 0, for all large enough N, 

card{sp(XAr) f]]p -e;p + e[} = kf, 

• for almost all rj > 0, for all small enough e > 0, for large enough N , 

sp(X^)flC\<C U ]p-e;p + e[. 

Remark 4.1. The proof of this theorem, starting at the end of this section, 
and completed in the next one, works without any change when the spikes may 
depend on N , more precisely if we assume that the spectrum of consists of 
N — r eigenvalues aj^' satisfying 

max distfa':^-', supp(/i)) — > 0, 

l<j<N-r ^ N^+oc 

and r spikes 6^ > ■ ■ ■ > 0^ such that: 

Vi e {fci + . . . + + 1, . . . , fci + . . . + fcj}, 6if — ^ Oj. 

The only reason we chose not to write the proofs under these slightly more 
general hypotheses is to avoid confusion in the notations. 

Remark 4.2. Actually, the conclusion of the preceding theorem holds for a 
random matrix Xj^ — Aj^ + , where A^ and Bn are independent random 
Hermitian matrices satisfying almost surely the assumptions given in the first 
section and in the preceding remark, under the extra assumption that the dis- 
tribution of the random matrix Bn is invariant by conjugation by any unitary 
matrix inV^. Note however that the quantities J,ki, . . . ,kj,9i, . . . ,9j are sup- 
posed deterministic. In particular, one recovers results from Ulf . U5f and U6] /. 
It is clear that a random Hermitian matrix Bj^ whose distribution is invari- 
ant by conjugation by any unitary matrix in has the same distribution as 
U'^BnUn, where B^ is a random real diagonal matrix, and Um is a random 
unitary matrix distributed according to Haar measure and independent of Bjq. 
The proof would then proceed, on the almost sure event on which A^ and B^ 
satisfy all the assumptions required, by conditioning with respect to the sigma- 
field generated by the sequences of random matrices A^, B^, and then applying 
Theorem \4.1\ 
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Example 4.1. When the rank of Apj remains finite, one recovers, using the 
preceding remark, results on the extremal eigenvalues from recalled in The- 
orem [n\ (see also the Gaussian case in \15\/ }. Indeed, in this setting, ii = Sq, 
and the set O is therefore the set of solutions p G R \ supp(i/) of equations 

as explained in Remark \3.S[ Theorem implies the conclusion of Theorem 
the condition 

{resp. 7r-i < 1/ lim G,y(z)) 

z-fa 

being equivalent to the existence of elements of O greater (resp. lower) than 
the maximum (resp. minimum) of the support of v, and the limiting points 
^^^(1/7*). (resp. G~^{l/-fr-])) being solutions of F^{p) = ji, (resp. F^{p) = 
"fr-j)- Notice that Theorem \4. 1\ deals in addition with the outliers of located 
in bounded components of the complement of the support of v. 

Example 4.2. When B]\j is drawn from the GUE(N ,a'^ ), which satisfies the 
assumptions stated in the preceding remark (invariance by unitary conjugation, 
almost sure weak convergence of the spectral measure towards the semicircular 
distribution J^, almost sure convergence to of the distance between the eigen- 
values and the semicircular support one recovers a particular case of the 
results on the outliers from '161. Indeed, in this setting, v is the semicircular 
distribution, which is 'Si -infinitely divisible, and the set O is therefore described 
by Remark \3.Sl The conclusion of Theorem \4.1\ is then exactly the one of the 
main result of U6f . 

Analogously, when Bj^ is a N x N Wishart matrix, which also satisfies the 
assumptions stated in the preceding remark (invariance by unitary conjugation, 
almost sure weak convergence of the spectral measure towards the Marchenko- 
Pastur distribution, almost sure convergence to of the distance between the 
eigenvalues and the Marchenko-Pastur support [Sfj ), one recovers the result 
on the outliers established by the two last authors of this paper and presented 
in I^OJ /. In this setting, v is the Marchenko-Pastur distribution, which is ffl- 
infinitely divisible, and the set O is therefore described by Remark \3.Si The 
conclusion of Theorem \4.1\ is then exactly the one of the main result of chapter 
7 oflMi- 

4.2 Reduction of the problem to the almost sure conver- 
gence of a r X r matrix 

In this section, we explain how we reduce the problem of locating outliers of 
Xm to a problem of convergence of a certain r x r matrix in the spirit of . 
Due to the invariance of the Haar measure under multiplication by any unitary 
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matrix, we may assume without loss of generality that both An and are 
real diagonal matrices: 




Diag(/3}^ 



Moreover, from the beginning of our argument, we will make use of the following 
additive decomposition of Ajy: 

An ^A'n + A'I^, 
A' TV / (^) (^) ^ 

A'^ = Diag(6ii - a, . . . , 6ii - a, . . . , 6*7 - a, . . . , 6lj - a, 0, . . . , 0), 



where the choice of a g supp(/i) is made so that limj^^o Gf^{iy+a) G M+i[— cxd, 0). 
Note that A'^ — *P9P, where P is the r x N matrix defined by 

P = {Ir\Orx{N-r))i 

Q is the r X r matrix 

9 = Diag(0i — a, ... ,01 — a, ... ,9 J — a, ... ,9 J — a), 



and *X denotes the transpose of the matrix X. Under our assumptions, the 
spectral measure of A'j^ (resp. Bn) weakly converges to (resp. ly), and all 
eigenvalues of A'j^ (resp. Bn) belong to any given neighborhood of supp(/i) 
(resp. supp(z/)) for large enough N gN. Hence, applying Corollary 3.1 of [15], 
one gets that, for any fc e N*, almost surely, 

3Nk e N,ViV > Nk,sp{A'N + U*nBnUn) C Kf. 

k 

We obtain thus the almost sure existence of a sequence {r]N)N>Ni of positive 
numbers converging to so that: 

ViV > A^i,sp(A^ + U*nBnUn) C K^^. 

(Choose for instance rjN = -^i for < N < TV^+i.) In the following, we will 
restrict ourselves to the almost sure event on which the sequence i'r]N)N>Ni is 
well-defined. Then, for iV > iVi, for any A G C \ K^^, the matrix A/jv - (A^ + 
UnBnUn) is invertible and 

dct(A/jv - Xn) = dct(A/jv - {A'n + U*nBnUn)) det{lN - i?Ar(A)*PeP), 
where 

RN{X) = iXlN-{A'N + U*NBNUN)r^ ■ (4.2) 
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Using that, for rectangular matrices X G Mjv,r(C),y € Mr^N{C), one has 
det(/Ar — XY) = det(/r — YX), one obtains: 

det(A/jv - Xn)) = det(A/jv - {A'j^ + Wj^BnUn)) det(/^ - PRN{\fPQ). 

Hence, for N > Ni, the eigcnvahies of X^ outside K^^^ are precisely the zeros 
of det(/r — Pi?jv(A)*P0) in that open set. We will denote by 

Mn := It - PRn^PQ (4.3) 

the analytic function defined on C \ K^^ , with values in the set of r x r complex 
matrices. 

Then, the following hindamental lemma allows to reduce the problem to the 
convergence of the sequence of analytic functions (Mjv)Ar>jVi- 

Lemma 4.1. Let M: C\K Mr{C) be a normal-operator-valued analytic 
function (i.e. M{z) G Mr{C) is normal for each z gC\K) so that 

(a) yzeC\K, M{z)* = M{z). 

(b) Qz>0 QM{z) invertible. 

Assume that there exists a sequence of positive numhers {riiq}^!^^ decreasing to 
zero and a sequence of analytic maps Mn: C \ K^^ — ^ Mr{C) so that 

1. there exists C > such that for all z G C such that \z\ > C, for any N, 

M]sf{z) is invertible, 

2. for any z € C\R we have Mn{z) G GLr{C). 

3. for any rj > 0, Mn converges to M , uniformly on C \ K^; 

If, for a fixed rj > such that the boundary points of are not zeroes of 
det(M), {pi, . . . ,Pp(n)} is the set of points z S C \ such that M{z) is not 
invertible, then 

(i) {pi,...,Pp(,,)} CM; 

(ii) dim(Ker(M(pj)) equals the order of zero of det{M {pj)) ; 

(Hi) For any < e < imin{|pi — pj\, d{pi,K^): I < i ^ j < there 
exists an iVo € N so that for any N > Nq the function det(MAr) is defined 
on C \ K^'', has exactly dim(Ker(M(pj))) zeroes in [pj — e, pj + e) for any 
j e {1, . . . ,Pp(r])} dfid exaciZj/ dim(Ker(M(pi)))+- • •+dim(Ker(M(/9p(^)))) 
zeros m C \ K^, counted with multiplicity, so that 

p{v) 

{zeC\K^: det(M;v(^)) = 0} C [j {pj - s, pj + e). 
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Proof: To begin with, by [H], M{z)M{z') = M{z')M{z)Jot all z,z' eC\K. 
Thus, there exists a unitary matrix U so that for all z G C \ K 

U*M{z)U = diag(/ii(z), . . . , hr{z)). 

Pick an X e^\K so that det M [x) = 0. The condition M{z)* = M{z) implies 
that M{x)* — M {x) whenever x € M, and thus hi{x) G R for any i e {1, . . . , r}. 
Let / ~ {ii, . . . , ij} be such that hi{x) = ii I £ I and hi{x) ^ else. Denote 
by TTifc the multiplicity of the zero of hi^{z) at x. The zero of det(M(a;)) is of 
order equal to mi + m2 + • • • + rrij. We only need now to argue that = 1 
for all 1 < fc < j. According to (b), ^hj{z) 7^ whenever Qz 7^ 0. This, in 
particular, implies that Qz'^hj{z) has constant sign on half-planes, and by the 
Julia-Caratheodory Theorem, h'j{x) 7^ for any x £ 'R\K. This proves (i) and 
(ii)- 

If we denote by f{z) — det(Af(z)) and fniz) = det(MAr(z)), then Hurwitz's 
Theorem jSTJ Kapitel 8.5] guarantees that for N large enough, /at will have 
exactly as many zeros - multiplicity included - as / has in C \ K'^ - and since all 
zeros of /at are known to be real, in C \ - and these zeros will cluster towards 
{pi, . . . , Pp(r{)\ in the sense that for any given e > there exists an e N so 
that 

{zeC\if«: det(Mjv(z))=0}c \^B{p^,e) 

i=i 

whenever N > N^. Moreover, for e > small enough, there are exactly 
dim(Ker(Af(/3j)) zeros of /^r in B{pj,s), multiplicity included. Since by 2., 
M^r is invertible in the two half-planes, we must have 

{zeC\K^: det(Mw(z)) = 0} C |J (p, - e, p, + e). 

□ 

5 Convergence of defined by ( 14.31) 

5.1 Preliminary results on the resolvent 

We begin this section by recording some facts on the resolvent Rn defined by 
(|4.2p . Recall that, if X is a selfadjoint operator on a Hilbert space with spectrum 
(j(X), then we shall denote by Rx{z) — [z — X]^^ its resolvent. It is known that 
this resolvent is analytic on C \ <y{X). If 1^9 is a positive unital linear functional 
on the unital algebra generated by X, then the distribution /ix,<^ of X with 
respect to can be recovered as 

G^,^^{z)^^{Rx{z)), zia{X). 
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Since in most cases it will be clear from the context which functional ip is 
considered, we shall suppress from the notation ^Jix,^p■ 

A more general notion of resolvent of X, which we shall use only sparingly 
in this paper, can be defined, following Voiculescu, as below: for an arbitrary 
operator h on the same Hilbert space as X, we can write its decomposition 




where 3?&, 96 are selfadjoint. We shall write 56 > if 36 > as operator 
on Hilbert space, and (36)^^ exists and is bounded. It has been noted by 
Voiculescu [S^ that 

Rx{h) = ih - X)-^ , 36 >0 (5.1) 

is an analytic map so that 3i?x(6) < 0. Moreover, as noted in [8, Remark 
2.5], if E is a positive unit-preserving linear map which leaves the algebra of 6 
invariant, then 

'^[E[Rx{b)]]-^ >^h. (5.2) 

For ah z e C \ R, Rn{z) defined by (g^ satisfies: 

VzeC\M,||i?jv(z)|| < (5.3) 

This implies that, for any z G C \ M, the matrix E[i?Ar(z)] is well-defined and 
also satisfies: 

||E[i?^(^)]|| < 

Lemma 5.1. When ^ 0, E[i?Ar(z)] is an invertible matrix and 

AC 

VzeC\M, ||E[i?w(z)]-i|| < |z| + Ci + — ^, (5.4) 

|!yz| 

where C I is any constant greater than sxvpj^f{\\A^\\ + \\BN\\) andC2 any constant 
greater than sup^ (trAr(_B^) — [ir jsi {B fq)]^^ . 

Proof. As noted in equation (|5.ip above applied to 6 = z — A'j^j, 5(z — A'j^ — 
U'^BnUn)~^ < 0. Since E is both positive and faithful, it follows that for 
any z £ C+, [(z — A'j^ — U^BnUn)^^] < 0, and thus by the same remark 
of Voiculescu Ji5|, E [(z — A'j^ — U^BnUn)"^] is invertible. The second state- 
ment of the lemma is equivalent to a statement about the power series expansion 
of z n- E[i?Ar(z)]^^ around infinity. The power series expansion 

E[i?.(z)].f:M^-±^l^^^, izi>p,ii+iii?..ii, 

n=0 
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assures us that 

f{z) {z~Aj^-E[U*^BNUN])-nRNiz)]-^ 



z-]K[A'r^ + U*^BNUN] ~ 



- 1 



E 



yfl + l 



1 

2; 

X 

1 r 

z 



E[(A^ + U*^BNUNf] - E[Aj^ + U*^BnUn]^ + - ■ 0(1) 



1 + E 



E[{A'j^ + U*^BNUNr] 



E[{Aj^ + W^BnUn?] - E[Af, + W^BnUn^ 



rO(l). 



Bwll is 



This power series expansion holds uniformly in N as long as ||A 
bounded uniformly in N . In particular, we obtain 

lim zf{z) = E[{Aj,+U*^BNUNf]-E[Aj,+U*^BNUN]'' ^trN{{B~tvN{B)f)-l, 

z— ^00 

uniform limit in N . 

As noted in equation (|5.2I) above, Qf{z) < 0, so for any positive linear 
functional (p on A'In{C), the function z H> (p{f{z)) maps C"*" into the lower half- 
plane, and \im^^^ zip{f{z)) = (p{l){tiN{B'^) - [tTN(B)]'^). Thus, z ^ (p{f{z)) 
is the Cauchy-Stieltjes transform of a positive measure supported on [—11^^^ 1 1 ~ 
IIBatII, P^ll + \\Bn\\] of total mass ip{l){tTN{B'^) - [tiNiB)]'^). It follows that 

\^{f{z))\ < ^{iTN{B^) - [tTNiBT), z e C+. 

Now, since positive linear functionals on von Neumann algebras reach their 
norm on the unit, the Jordan decomposition of linear functionals allows us to 
write 

||/(z)|l < sup |^(/(z))| < 4 sup |^(/(z))| < ^{tT^iB') - [tTNiB)n 

llvll = l ¥'>0,v(l) = l 

Since E[Rn{z)]^^ = z — E[A^ + [/]^_Bjvt^Jv] — /(-z), for any z e C\M, we conclude 
that 



\\E[Rn{z)]-^\\ < \z\ + \\Aj,\\ + \\Bn\\ + \\f{z)\\ < \z\ + Ci + 
as stated in our lemma. 



□ 
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One now states concentration results that will allow to reduce the proof of 
Proposition 15.11 to the convergence of a sequence of deterministic matrices and 
to estimate the variance of each entry of the resolvant Rn{z). 

Lemma 5.2. (i) Vz e C\M,Pi?jv(^)*-P- ™[i?Ar(z)]*P ^ 0. 
(ii) Vz e C\R.y{k,l) e {1, . . . ,N}\Y{{Rn{z))m) < 

Proof. Fix z S C \ M. It is clear from the remark that writing P{Rpf{z) — 
E[Rn{z)]YP corresponds to taking the upper left r x r corner of Rn{z) — 
M[Rn{z)], that (i) is equivalent to: 

V(fc, e {1, . . . , rf, {Rn{z) - E[i?^(z)]),, 0. (5.5) 

Now, for any (fc, I) € {1, . . . , iV}^, since the function 

f -.Un ^ RN{z)ki 

is Lipschitz on the unitary group Uat with Lipschitz bound j^p, by Corollary 
4.4.28 of the book [1 , for any Q < a < \, 

P {\ {Rn{z) ~ mN{z)])^i I > < 2exp (-CiV2"|5z|4e2) . 

Hence, one gets (|5.5p by a standard application of Borel-Cantelli lemma, and 
{a) by the classical formula holding for a positive random variable X: 



E(X) = / P(X > t)dt. □ 
"'0 



5.2 Convergence of Mjv 

We investigate the convergence of the sequence of analytic functions {Mm)n>Ni 
defined by (|4.3p . We need the following preliminary lemma. 

Lemma 5.3. The function x : z i— > is analytic on C \ supp(/i ffl ly) and 

satisfies xiz) — x{z) for any z in C\ supp(/i ffl v). 

Proof. This lemma readily follows from Lemma |3. II □ 

Proposition 5.1. Almost surely, for any rj > 0, the sequence {Mn)^^^^, where 
No > -^1 is such that ViV > iVo, rjpj < rj, converges to 

M = diag(l - (01 - a)x, . . . , 1 - (0j - a)x), 

uniformly on compact subsets of C \ K^. 
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Since 

Mat ^Ir- PRn PQ, 
it is equivalent to prove the convergence of PRn^P towards x^r- 

We first study the convergence of the sequence of analytic functions (E[i?Ar])7v>i : 
Proposition 5.2. 

VzeC\M, PE[Rn{z)YP — > x{z)Ir- (5.6) 

Fix z G (which is sufficient by a reflection argument). We now break the 
proof in three lemmas. First, a strenghtening of the result from [22] stating that 
that the matrix ]E[i?Ar(z)] is diagonal: 

Lemma 5.4. // 6 £ M]y(C) is so that b — U^BnUn is invertible for each 
value of the Haar unitary Un, then E [(6 — U'^BnUn)^^\ G where the 

hicommutant {5}" is taken in Mjv(C). 

Proof. Pick an arbitrary unitary V in the commutant {6}' of b. By the invari- 
ance of the Haar measure, 

V*E[{b-U'^BNUNy^]V = ¥.[{V*bV -V*U*i^BnUnVY^] 

= E[{b-U*j,BNUN)-^], 

so that E\{b — UIjBis!Un)~^\ G {V}' ■ Thus, since a von Neumann algebra 
equals the span of its unitaries, E [(6 — U^BmUn)~^\ € {b}", as claimed. □ 

Recall that the first r eigenvalues of A'j^ are all equal to a. We apply the above 
lemma to 6 = z — A'j^ to conclude that the first r eigenvalues of E [i?jv(z)] are all 
equal, and thus PE[Rpf{z)YP — XN{z)Ir for the Cauchy-Stieltjes transform 
Xn of some probability measure depending on A'^, B^y. 

Our next task is to establish an approximate matricial subordination equation, 
namely to prove that W.[Ri^{z)] is asymptotically equal to {ujn{z)In — ^jv)~^i 
for a certain complex number ujn{z). Then, we prove the uniform convergence 
on the compact subsets of C+ of the sequence of analytic functions {ujn)n>i 
towards u. 

Lemma 5.5. For z G C+, one has: 

\\E[RN{z)]-{LUN{z)lN-A'r^)-^\\ 0, 

A/"— f+oo 

where ^ ^ 

E[i?Ar(z)Jii Xn{z) 
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Proof. First notice, using (|5.2p . that ujm defined by (15.71) . satisfies: 

Vz e C+, 3wjv(z) > 32. (5.8) 

Fix z £ C+ and define 

^Niz) :=E[i?jv(z)]-i+A^, 

which belongs, according to Lemma to {A'^}". We denote its fc-th diagonal 
entry 

E[i?Ar(z)Jfefe 

Note that ujn{z) — {VLn{z))xi and that, using (|5.3p and ([5 



and 

If we prove that: 

3C'{z) > 0, ViV > 1, V(fc, /) e {1, . . . ,iV}, - (r!Ar(z))z,| < 

then we may conclude: 

\\¥.[Rn{z)]~ {ujn(z)In - A'j,)-^\ 

= m^iz) - - {lon{z)In - A'nVW 

\mN{z) - A'j,)-\ujn{z)In - nN{z)){uJN{z)lN - A'j,)-^\\ 
< \mM{z) - A'j,)'^\\\\{ojn{z)In ~ nN{z))\\\\{u^N{z)lN - A'j.r^W 

^ T^\mN(z)UlN~n^{z)\\ 



< 



|^Z| 

C'{z) 





For (fc, /) e {1, . . . , A^}, observe that 

{SlN{z))kk - (r!iv(^))« = {^n{z)Eu - EkinN{z))ki. (5.9) 

Define, for a given deterministic matrix X e Mn{C), 

Ax := E[{Rn{z)-E[Rn{z)]){A'j^X -XA'^)iRN{z)~E[RN{z)])] 

= E[Rn(z){A'j^X - XA^)i?Ar(z)] - E[Rn{z)]{A'j^X - XA'j^)E[Rn{z)] 
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Noting that for any deterministic Hermitian matrix X, 

jE {{z -A'j,- e-**^'C/^i?w[/we"^)-i)|^^^ = 0, 

we readily deduce that 

E[Rn{z){A'j^X - XA'j^)Rn{z)] = E[Rn{z)]X - XE[Rn{z)1 

and then extend this identity by linearity to any matrix X G Afjv(C). It follows 
that 

nNiz)X - X^Niz) = -E[Rn{z)]-^AxE[Rn{z)]-\ (5.10) 
For X ^Eku and (fSAO]) yield 

\(SlN{z))kk - [^n{z))ii\ 

= \{E[Rn{z)]-^I^e,MRn{zT%i\ 

< \E[RN{z)]^^\\{AEjkmRN{z)]-'\ 

X E[\{Rn{z) - E[RN{z)])kk{RN{z) - E[Rn{z)])u\] 

< 2\\A'j,\\(\z\ + C, + ^) Y{iRMiz))kk)^Y{{RNiz))u)-- 



< 



Ci\z\ + l+^ 
N\^z\* 



,-JZ\ 

2 



where we used (|5.4[) and Lemma [5.21 (ii) in the three last inequalities. And we 
are done. □ 

We now study the convergence of the sequence {ujn)n>i defined by (|5.7p . 

Lemma 5.6. The sequence of analytic functions {ujn)n>i defined on C \ R 
converges uniformly towards lu on the compact subsets o/C \ R. 



Proof. It follows from Lemma 15.51 by taking the normalized trace, and using 
the notation 

QNiz) := E\Ga , . (z)] 

that 

gN{z)-Gf, {u:n{z)) 0. 

N N^ + OO 

Using Lemma 7.7 in [T3] and (|5.3p . we deduce that 

gN{z)-G,{uJN{z)) 0. (5.11) 

N^-\-oo 
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The sequence of analytic functions (wjv)Ar>i is normal, and thus there exists 
at least one converging subsequence. For any fixed z S C+, let us consider a 
converging subsequence a;0(jv)(z) of L>Jf^{z) and denote by l{z) the limit. As 
noted above in (j5.8p . ^ixJt^{z) > 5z. Thus, if '^z is large enough, |a;Ar(z)| will 
be large, and we can then uniquely invert with respect to composition uni- 
formly in N in (|5.11l) to obtain G~^{gN{z) — o(l)) — ujn{z). Letting TV go 
to infinity, Voiculescu's asymptotic freeness result guarantees that gN{z) — > 
Thus, l{z) = G'^^ {G ^^y{z)) , independent of the convergent subse- 
quence we chose. This implies that lim7v->.oo = ' uniformly on 
compact sets of C+, and by analytic continuation that I = oj. Therefore, for 
any z £ C"*", uj{z) is the unique cluster point of itj]\[{z). □ 

Proof of Proposition 15. 2i Fix z E C+ (which is sufficient by a reflection 
argument). We proved that 

PE[Rn{z)YP = -^^ Ir. (5.12) 

We simply conclude by using Lemma 15.61 and '^ujn{z) > 3z, 3a;(z) > 3z. □ 

It will be important in the proof of Proposition l5.1l to have the following analytic 
continuation result. 

Theorem 5.1. (lemma 7.6.5 \17l ) Let A he an open nonempty subset of C 
and D <Z A such that D = A. Let {gn) he a sequence of locally hounded holo- 
morphic functions on A such that lim„_j.+oo gn{d) exists for any d E D. Then, 
the sequence (gn) converges towards a function g which is holomorphic on A, 
uniformly on each compact suhset of A. 

Proof of Proposition 15. 1[ Define 

■D,^ = {zeC\ K^, nz e Q, 3z e Q*}. 

According to Lemma [5?2] (i) and (|5.6p . for any z G Vr,, almost surely, PRn{zYP 
converges towards xiz)Ir. Hence, almost surely, for any fc, Z G {l,...,r}, for 
any ry > 0, {RN)ki is a bounded sequence of holomorphic functions on C \ 
such that the limit of {RN{z))ki exists for any z G Vrj. Therefore, according to 
Theorem 15. 11 we can deduce that, almost surely, {RN{z))ki converges towards 
an holomorphic function Xki on C \ K'^, uniformly on each compact subset of 
C \ K^. Of course, Xki coincides with ^^^^^^ on so that Xki — &kiX- The 
proof is complete. □ 

We will now prove Theorem 14.11 by applying Lemma 14.11 on an almost sure 
event on which its assumptions are satisfied. 

Proof of Theorem 14. 11 We consider the almost sure event, whose existence is 
guaranteed by Proposition 15. 11 on which there exists an integer A'^i, a sequence 
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{i1n)n>Ni of positive numbers converging to 0, so that 

and, for any 77 > 0, the sequence {Mn)]^^]^^, where Nq > Ni is such that 
V-/V > 7 Vn < converges to 

M = diag(l - (01 - «)x, . . . , 1 - {ej ~ a)x), 

uniformly on compact subsets of C \ K^. On this event, apply Lemma 141] to 
the sequence {M]s[)n>No and its uniform limit M. The funtion M: C \ K ^ 
Mr{C) is indeed a normal-operator- valued analytic function satisfying trivially 
conditions (a) and (b) of Lemma |4?T] The sequence {Mn)n>No consists of 
(random) analytic maps on C\ K^^. Condition 3. of Lemma 14. II is guaranteed 
by Proposition [23] Condition 2. is straightforward. To check condition 1., it is 
sufRcient to argue that for any z such that \z\ > Ci, 

\\Rn{z)\\ < 



d(z,hCi,Ci])' 

where Ci is chosen as in Lemma 15.11 Almost every 77 > is such that the 
boundary points of are not zeroes of det(M), so for such 77's, we get exactly 
the conclusion of Theorem 14. II Indeed, as explained in Section 3.2, eigenvalues 
of Xn in C \ are exactly zeroes of det(MAr), and the set of points z such 
that M{z) is not invertible is precisely O. □ 
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